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Abstract 
Ground-state and thermodynamic properties of the one-dimensional 
^ Heisenberg antiferromagnet in which two S* = 1/2 and two S = 1 spins 
are arranged alternatively are studied by a quantum Monte Carlo method 
and by analytical estimates. It is found that the ground state of the system 
undergoes a second-order phase transition which accompanies the vanishing 
of the energy gap. 
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One-dimensional quantum Heisenberg antiferromagnets have been the 
subject of numerous studies over a long period, in particular since Haldane's 
prediction [1,2] of the difference between the integer-spin and half-integer- 
spin cases. Almost all results of these studies support the prediction, and 
now it is widely agreed that for isotropic nearest-neighbor (nn) exchange, 
the former case is massive, while the latter case is massless. Under these cir- 
cumstances, the 5*= 1/2 and S=l alternating spin chain with isotropic nn 
exchange, which is a simplest example of mixed spin systems, has recently 
been of considerable interest both experimentally [3] and theoretically [3-7]. 
A characteristic feature of this chain is that it is a quantum ferrimagnet 
and its ground state is magnetic. 

In order to explore the mixed spin chain with a nonmagnetic ground 
state, we discuss here the one-dimensional antiferromagnet in which two 
S* = 1/2 and two S = 1 spins are arranged alternatively. We express the 
Hamiltonian describing the system as 

N/4 

'H = ^('/l S4^-3-S4^-2 + -^2 S4^-2-5'4£_i + J3 Su-l-Su^J2 Su-SU-\-l) , (1) 
i=l 

where si {i=l, 2, 5, 6, • • •, A^-3, N-2) and Si (f = 3, 4, 7, 8, • • •, N-1, 
N) are, respectively, the 5=1/2 and S = l spin operators, and N, being 
assumed to be a multiple of four, is the total number of spins. We impose 
periodic boundary conditions (sjv+i =si). It is noted that the ground state 
of the system is nonmagnetic except when both Ji and J3 are negative. 

Applying the Wigner-Eckart theorem, we may conclude that in the 
limit of J2 ^ 00 ( J2 -^ —00), the present system is equivalent to the 5* = 1/2 
{S = 3/2) bond-alternating chain with the exchange constants | Ji and 
^-^3 (^-^1 and I J3) as far as sufficiently low energy states are concerned. 



Thus, the system is gapless when J2^oo and Ji = 16J3>0 and also when 
J2 ^ oo and Ji = 4 J3 > 0. A perturbation calculation taking into account 
the lowest singlet and triplet of each unit cell (54^-2, S^i-i^ Su, 54^+1) 
leads to the result that, for sufficiently small values of |Ji| and | J2I with J3 
fixed at 1.0, the energy gap vanishes when Ji = | J2^ + 2J2^ + - • •. 

Confining ourselves to the case where Ji = J3 = 1.0 and J2>0, we in- 
vestigated numerically the ground-state and thermodynamic properties of 
the system by a quantum Monte Carlo (QMC) method. Let us first discuss 
the ground-state properties. We carried out, at sufficiently low tempera- 
tures (T = 0.05 and 0.10), a QMC calculation without the global fiip within 
the subspaces of M = ^^^( {s\^_^^ s\^_2^ Sl^_-^^ Sl^ =0 and 1 to esti- 
mate the energies E[M) for M = and 1, the energy gap A = £'(l)— £"(0), 
and the nn two-spin correlation functions ui^2 = (sfsl), 0J2,3 = (s^'S'l), and 
a;3^4 = {S^Sl), where (• • •) stands for the Monte Carlo average within the 
M = subspace. In the calculation, 10^ Monte Carlo steps were performed 
after 10^ initial steps for obtaining the thermal equilibrium. The Trotter 
number n was chosen to be n = 12, 16, 24, 32, 40, and 48, and the n- 
dependence of the QMC results was extrapolated to n ^ 00 by making a 
least-squares fit to the formula, An = Aoo^ai/'n?^a2/n'^- For both temper- 
atures of r = 0.05 and T = 0.10, almost the same extrapolated values were 
obtained, and thus we consider that T = 0.05 is low enough to discuss the 
ground-state properties. The QMC calculation was carried out for A^ = 8, 
16, 32, 64, and 128, for various values of J2. As for £"(0), £"(1), and A for 
N = 8, we compared the extrapolated values for T = 0.05 with the values 
obtained by a method of exact diagonalization; this comparison shows that 
the precision in the former values goes down to the second decimal point. 
In order to obtain the results in the thermodynamic limit {N -^ 00), the 



extrapolated values for the latter four N's were further extrapolated by the 
least-squares method using a linear function of l/N"^ . 

Fig. 1 depicts the J2-dependence of the energy gap A. From this fig- 
ure we see that there exists a massless point at J2 = J2c~0.77, at which 
the ground state of the system undergoes a second-order phase transition. 
Giving a thought to the J2 -dependences of the correlation functions a;i,2, 
0^2,3 , and 0^3 ^4 depicted in Fig. 2, we may schematically represent, by means 
of the valence-bond-solid (VBS) picture [8,9], the ground states for J2 < J2C 
and that for J2 > J2C as shown in Fig. 3. 

Next, we briefly discuss the thermodynamic properties. Carrying out 
a QMC calculation which takes the global flip into account, we calculated 
the temperature dependence of the magnetic susceptibility and the spe- 
cific heat. We obtained the result that the susceptibility per unit cell for 
J2 = 0.77 takes a finite value in the T ^ limit; this is consistent with 
the fact that the energy gap vanishes at this value of J2. We also found 
that for sufficiently large values of J2 such as J2 = 5.0 and 10.0, the spe- 
cific heat per unit cell versus temperature curve has a double peak; the 
higher-temperature peak is associated with the Schottky-type peak of the 
two-spin system oi S = 1/2 and S =\ spins, while the lower-temperature 
one is associated with the specific heat peak of an antiferromagnetic 5* = 1/2 
chain with bond alternation. 
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Figure Captions 

Fig. 1. Plots of the N^oo energy gap A versus J2. The open circles show 
the QMC results, and the solid lines are guides to the eye. Note that 
A vanishes at J2 = J2c ~ 0.77. 

Fig. 2. Plots of the N ^ 00 correlation functions 0^1,2 (open circles), 0^2,3 
(open squares), and 0^3,4 (open diamonds), versus J2, where each func- 
tion is normalized by its minimum value (note that — l/4<a;i^2 ^ 1/12, 
— 1/3 < 0^2,3 < 1/6, and —2/3 < us,4 < 1/3). The solid lines are guides 
to the eye. 

Fig. 3. Schematic representations of the ground states for (a) J2 < J2C and 
(b) J2 > J2C by means of the VBS picture. The solid circles represent 
the 5* = 1/2 spins, and two 5*= 1/2 spins connected by the solid line form 
a singlet pair. Each open ellipse surrounding two S = 1/2 represents 
an operation of constructing an S' = l spin from these 5* = 1/2 spins by 
symmetrizing them. 
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